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18 Abstract 

19 Standing genetic variation and the historical environment in which that variation arises 

20 (evolutionary history) are both potentially significant determinants of a population's 

21 capacity for evolutionary response to a changing environment. We evaluated the relative 

22 importance of these two factors in influencing the evolutionary trajectories in the face of 

23 sudden environmental change. We used the open-ended digital evolution software Avida 
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24 to examine how historic exposure to predation pressures, different levels of genetic 

25 variation, and combinations of the two, impact anti-predator strategies and competitive 

26 abilities evolved in the face of threats from new, invasive, predator populations. We show 

27 that while standing genetic variation plays some role in determining evolutionary 

28 responses, evolutionary history has the greater influence on a population's capacity to 

29 evolve effective anti-predator traits. This adaptability likely reflects the relative ease of 

30 repurposing existing, relevant genes and traits, and the broader potential value of the 

31 generation and maintenance of adaptively flexible traits in evolving populations. 
32 

33 Introduction 

34 The diversity and complexity of any biological system reflects past evolutionary 

35 responses to environmental conditions, as much as contemporary ones. Underlying those 

36 responses are a number of extrinsic and intrinsic factors influencing populations' capacity 

37 to evolve responses through the generation or utilization of genetic variation [1,2]. 

38 Because evolutionary responses build on available heritable variation, mechanisms that 

39 influence the generation and maintenance of variation [3] can strongly shape the 

40 evolutionary trajectories of populations [4]. Among the potential factors involved, 

41 standing genetic variation (SGV) and evolutionary history (EH) are likely to be 

42 significant determinants of adaptability to novel environments [5]. Accordingly, 

43 understanding how these critical factors either promote or constrain population 

44 evolutionary potential provides insight into the realized pathways that led to historic 

45 evolutionary outcomes, as well as those that will shape future populations. 
46 
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47 Standing genetic variation is the presence of alternate forms of a gene (alleles) at a given 

48 locus [5] in a population. While an allele may be mildly deleterious or confer no fitness 

49 advantage over other forms under a given set of environmental conditions [6], that allele 

50 may become beneficial if the environment changes. As selection can act only on available 

51 variation, SGV provides a potential means for more rapid adaptive evolution [7] 

52 compared with the de novo mutations [5,8], particularly if environmental conditions 

53 change. 
54 

55 In addition to SGV, a population's historical selection environment (i.e. evolutionary 

56 history) may play a strong role in determining the speed and the extent to which 

57 populations can adapt to directional environmental change [9,10]. In particular, EH will 

58 have influenced the genetic variation and genetic architecture of traits in contemporary 

59 populations. If changes in the environment alter the strength of selection on a trait, 

60 populations with an evolutionary history of adaptation to similar pressures may be 

61 mutationally "closer" to the discovery of new [1 1], or rediscovery of historic beneficial 

62 traits [12]. 
63 

64 With the contemporary rise of experimental evolution as a means of testing evolutionary 

65 and ecological hypotheses [13,14], a valuable and untapped opportunity now exists to 

66 elucidate the roles SGV and EH have played in determining historic, realized rates of 

67 adaptive evolution. Furthermore, while SGV is known to be an important determinant of 

68 the speed of evolution [5], it is less clear how levels of SGV and EH, alone or in 

69 combination, impact the overall evolutionary potential of populations. Such an 
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70 understanding could allow population geneticists to evaluate how SGV and EH may 

71 contribute to or constrain the future evolvability of populations, particularly in human- 

72 modified environments [1, 15, 16] 
73 

74 To evaluate the individual and interactive effects of SGV and EH on evolutionary 

75 outcomes, we measured their relative and independent contributions toward evolutionary 

76 potential in populations using the digital evolution software Avida [17]. Digital evolution 

77 experiments carry several significant advantages for addressing evolutionary questions 

78 requiring systematic manipulation and highly controlled environments. Among these, 

79 generation times are rapid, the experimenter has full control over initial environmental 

80 conditions, and detailed genetic, demographic, and behavioral trait data can be recorded 

81 perfectly. Organisms in Avida can also engage in ecological competition and other 

82 complex interactions, and the system can allow for co-evolution with predators in 

83 predator-prey systems [18,19,20]. Furthermore, unlike in evolutionary models and 

84 evolutionary simulations that impose artificial selection via explicit selection functions, 

85 Avida uniquely allows for unrestricted, unsupervised, and non-deterministic evolution via 

86 natural selection [17]. 
87 

88 Predation is an ecologically important agent of selection [21,22] as demonstrated by the 

89 diverse array of prey defenses that have evolved [23,24,25,26]. Accordingly, we used 

90 Avida to test whether historic exposure to predation influenced how prey populations 

91 responded to pressures from new invasive predators. We then further examined which 

92 factors (SGV, EH, and their interaction) were important in determining the future 
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93 defensive and competitive abilities of those populations. While both factors played 

94 significant roles, we show that evolutionary history has the stronger influence on 

95 evolutionary potential, with the strength of SGV's role being contingent on a 

96 population's evolutionary history. 
97 

98 Materials and Methods 

99 Avida 

100 The Avida digital evolution platform is a tool for conducting evolutionary experiments on 

101 populations of self-replicating computer programs, termed digital organisms [17]. Digital 

102 organisms are composed of a set of instructions constituting their "genome". Organisms 

103 execute genome instructions in order to perform actions such as processing information, 

104 and interacting with their environment, and for reproduction. Additionally, predefined 

105 combinations of instructions allow organisms to consume resources from the simulated 

106 environment. Sufficient consumption of resources, to a level defined by the experimenter, 

107 is a prerequisite for organisms to copy their genomes and divide (i.e. reproduce). During 

108 the copy process, there is an experimenter-defined probability that genomic instructions 

109 will be replaced with a different instruction (substitution mutation) randomly selected 

110 from the full set of all 60 possible instructions (instruction set files are in the skel 

111 directory in the GitHub repository [ https://github.com/fishjord/avida_predation_scripts ]). 

112 Separately, there are also set probabilities for the copy of an instruction to fail (creating a 

113 deletion mutation) or for a chance insertion of a new instruction into the copy's genome. 

114 In real-time, generation times in Avida are typically only a few seconds, and 

115 experimental environments are fully definable. 
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116 

117 Environment 

118 We ran all Avida trials in a 251 x 251 bounded grid-cell environment. Each cell started 

119 with one unit of resource. If a prey organism fed in a given cell (by executing an 'output' 

120 instruction), the resource in that cell was completely consumed and that cell was 

121 replenished with resources over 100 updates. 
122 

123 Timing 

124 All Avida trials were run for a set number of updates, as indicated in each of the Phase 

125 descriptions below. Updates are the measurement of time in Avida, reflecting the number 

126 of genomic instructions executed across the population. Here, each organism executed 30 

127 instructions per update. Organisms had a maximum lifespan of 15,000 instructions 

128 executed (500 updates). Realized average number of generations for the populations in 

129 each treatment is given in their respective sections below. 
130 

131 Reproduction 

132 Organisms, both predator and prey, were required to have consumed 10 resource units 

133 and be at least 100 updates old to successfully reproduce. Genome copying and 

134 reproduction occurred when the organism executed a single reproduction ('repro') 

135 instruction. Execution of the reproduction instruction took one update, compared to 

136 1/3 0th of an update for other instructions (except for predation handling time, as noted 

137 below). When an organism reproduced, its offspring was placed in the cell the parent was 

138 facing. This new organism's genome was a copy of the parent's genome, with the 
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139 exception any substitutions, insertions, and deletions that occurred during the copy 

140 process. Substitutions were independent and only one insertion or deletion was allowed 

141 per reproduction (at rates specified in the 'Phase 1 ' and 'Phase 2' descriptions, below). In 

142 addition to the birth of the new organism upon reproduction, the parent organism was 

143 effectively reborn, with all internal states reset. Unlike in legacy Avida experiments, 

144 multiple organisms could occupy the same cell [as in 18-20]. Thus, newborn organisms 

145 did not replace existing occupants of the cells into which they were born. 
146 

147 Predation 

148 In Avida, predation occurs via the evolution of an 'attack-prey' instruction, the execution 

149 of which allows an organism to attack and kill a non-predator in the cell it is 

150 facing. Predation in Avida and the base predator Avida configuration used is described 

151 in [17]. Here, as handling time, the attack-prey instruction costs the predator 10 cycles 

152 (one-third of an update) if the attack is successful. Predators receive 10% of each prey 

153 individual's consumed resources after a successful kill, which is then applied toward the 

154 consumption threshold required for reproduction. An organism is classified as a predator 

155 in Avida if it attacks a prey organism, not simply by evolving the 'attack-prey' 

156 instruction. 
157 

158 In simple (i.e. lacking topographic features like refuges) and confined environments like 

159 the one used here, Avida trials containing predators typically require a minimum prey 

160 population level below which predator attacks fail. In top-down limited populations like 

161 those used here, minimum prey levels prevent population extinctions and serve to 
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162 standardize prey population sizes and, thus intra-specific competitive pressures. In 

163 practice, because prey are constantly being born into populations, prey kills are prevented 

164 for very short periods of time and failures simply require predators to make multiple, 

165 rapid attacks. Here, attacks on prey were always fatal if there were more than 1,000 prey 

166 in the environment (i.e., 1,000 was a minimum prey level). Similarly, in trials without 

167 predators, prey population sizes were controlled via a preset maximum population cap 

168 (1,000 organisms), as described for each Phase below. When a birth caused that limit to 

169 be exceeded, a random prey organism was removed from the population. 
170 

171 Phase One: Evolutionary History 

172 Our study was divided into three experimental phases (Figure 1). In the first phase 

173 ("Phase 1"), we evolved two sets of 30 base populations. Each of these populations was 

174 initiated via the placement of 9 identical prey organisms into a new grid environment 

175 (Figure 1, Phase 1). These ancestors had simple genomes that allowed them only to 

176 blindly move around the environment, repeatedly attempting to feed and reproduce. In 

177 one set of populations, the attack-prey instruction was allowed to appear in the 

178 organisms' genomes through random mutations, thus allowing for the potential evolution 

179 of predators (Figure 1, Phase 1). Henceforth we refer to these populations as "predator 

180 EH" treatments. In the second set of populations, the "attack" instruction was prevented 

181 from mutating into organisms' genomes. This second set of "no predator EH " base 

182 populations thus evolved in the absence of predators in Phase 1. Each population 

183 evolved for two million updates (on average, 20,000 generations). During this phase, 

184 reproduction carried a 25% probability of a single, random, genomic instruction 



8 



Downloaded from http://biorxiv.org/on September 18, 2014 

185 substitution (maximum one substitution per genome), and 5% probabilities (each) of 

186 single insertion and deletion mutations. 
187 

188 Phase Two: Standing Genetic Variation 

189 The second phase (Figure 1, Phase 2) of the experiments allowed us to evaluate how the 

190 amount of SGV as well as EH influenced the evolutionary potential of prey in the 

191 presence of novel predators. For all Phase 2 trials, prey populations (compositions 

192 described below) were evolved for 200,000 updates in the presence of identical copies of 

193 a single static predator population (i.e. predators in Phase 2 were prevented from 

194 evolving further, dying, or reproducing). Taking in sensory information is a prerequisite 

195 for behaviorally responding to and targeting prey, and we used visual input rates as a 

196 proxy for predator effectiveness. If organisms evolve to use visual sensors, these sensors 

197 provide them with information about their environment, allowing them to "look" at the 

198 objects (food, prey or other organisms) in the 45-degree area in front of them. Because 

199 predator populations varied across Phase 1 replicates, we selected the Phase 2 predator 

200 population from the Phase 1 replicate with highest realized usage of visual sensors (i.e. 

201 highest per-predator 'look-ahead-intercept' instruction executions) to ensure that an 

202 effective predator population was used in Phase 2. Ultimately, the predator population 

203 selected and copied for use in all Phase 2 trials had 209 predator organisms. 
204 

205 To create Phase 2 Evolutionary History (EH) treatment source populations, we first 

206 selected a single prey population from a random Phase 1 replicate in each EH treatment 

207 {predator and no predator EH), excluding the population from which the Phase 2 
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208 predator population was drawn. For each of the two selected populations, we excluded 

209 any organisms that were classified as predators or had parents that were predators (as 

210 determined by an internal state, see [18,19]), to ensure that Phase 2 prey populations 

211 consisted only of prey organisms. An organism must execute an attack instruction to be 

212 designated a predator, some prey organisms may have attack instructions in their 

213 genomes, provided they have not yet used them; we therefore replaced any attack 

214 instructions in the remaining prey populations with a 'do nothing' (nop-X) instruction, 

215 and prevented any new copies of the attack instruction from mutating into offspring 

216 genomes. 
217 

218 To create Phase 2 SGV treatment populations, we created separate "clone ", 

219 "intermediate" , and "high " SGV populations by sampling each of the two EH treatment 

220 source populations ("seed populations"; one predator EH population and one no predator 

221 EH population: Figure 1, Phase 2). First, for high SGV populations, we simply used 

222 duplicates of the selected EH treatment source populations. Second, we created each 

223 clone SGV population by randomly selecting a single genotype from a seed population 

224 and making as many duplicates of that genotype as there were organisms in the seed 

225 population. A genotype could be selected for multiple Phase 2 clone replicates. Third, 

226 we created the intermediate SGV populations by randomly sampling genotypes (with 

227 replacement) from a seed population, creating up to 55 identical copies of each sampled 

228 genotype. We repeated this process until the new intermediate SGV population was the 

229 same size as the seed population; this sampling scheme yielded intermediate SGV prey 

230 populations the same size as their high SGV seed populations, but with -50% lower 
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231 diversity (Shannon's diversity index), on average. In all, there were six (2 EH x 3 SGV) 

232 Phase 2 treatments, each with 30 replicates. 
233 

234 For each of the six combinations of SGV {clone, intermediate, high) and EH {predator, 

235 no predator EH), we evolved 30 replicate populations for 200,000 updates (500 

236 generations, on average), each in the presence of a copy of the constructed Phase 2 

237 predator population. The mutation rates in Phase 2 were lowered to 0. 1% substitution 

238 probability, and 0.5% insertion and deletion probability. We found empirically that these 

239 mutation rates yielded, on average, 5% divergence (see below) between ancestor and 

240 final organisms in Avida after 500 generations (200,000 updates using the Phase 2 

241 settings), matching the expected divergence of a bacterial genome over 500 generations 

242 [27]. Other configuration settings were identical to those used in Phase 1. 
243 

244 At the end of Phase 2, we calculated how different the resulting population was from the 

245 starting population by using a dynamic programming algorithm [28]. We measured the 

246 genetic divergence of each prey organism in each population by aligning the current 

247 genome sequence with that of its ancestor from the beginning of Phase 2. Each organism 

248 from the seed population was tagged with a unique Lineage ID, which is shared with all 

249 progeny of the organism. The Lineage ID of the organisms present in the population at 

250 the end of Phase 2 could then be used to identify their ancestor from the beginning of 

251 Phase 2. We then calculated the divergence as the percent identity between the two 

252 aligned genome sequences. Alignment was required since genome sizes changed during 
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253 Phase 2, via insertions and deletions. 
254 

255 Phase 3: Competitive Evaluation of prey populations 

256 For the third and final phase (Figure 1, Phase 3a), we used a set of both "ecological" 

257 evaluation simulations to measure the fitness of the final prey populations in the presence 

258 and absence of predators, and trait assays to evaluate the evolution of anti-predator traits 

259 during Phase 2. For all Phase 3 evaluations, substitution, insertion and deletion mutation 

260 rates were set to zero. For all types of Phase 3 trials (described below), in order to create 

261 an uneven resource landscape (i.e. in the absence of consumption by prey, all cells would 

262 contain resources that could be consumed), populations were introduced into their test 

263 environments, run for an initial 1,000 updates, and then reintroduced a second time. Once 

264 prey were reintroduced, we ran each trial for an additional 10,000 updates, and recorded 

265 population sizes for the two competing populations. 
266 

267 As in Phase 2, Phase 3 trials used copies of a single Phase 1 predator population. The 

268 Phase 3 predator was intended to exhibit novel strategies relative to the Phase 2 predator 

269 population. We selected the Phase 3 predator population by first eliminating the three 

270 Phase 1 replicates used for creating Phase 2 predator and prey populations, and then 

271 selecting one predator population that exhibited an average level of visual sensor usage 

272 (about half that of the predator population selected for Phase 2). The final selected 

273 predator population for Phase 3 contained 253 predators. 
274 
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275 Phase 3 involved two sets of trials: the first set allowed for evaluating competitiveness, 

276 while the second set allowed for evaluating the effects of SGV and EH on expressed 

277 levels of evolved anti-predator behavioral traits. In the first set of Phase 3 trials (Figure 1, 

278 Phase 3a), we conducted pairwise competitions in which each population from each of 

279 the six Phase 2 treatments was paired, with a randomly selected population from each of 

280 the five other treatments (for a total of 30 pairings x 30 replicates = 900 competitions). 

281 For each pairing, the two populations were injected once into an environment with the 

282 Phase 3 predator population ^predators present PT ") and once into an environment 

283 without predators ( "predators absent PT "). For each of these trials, we enforced a prey 

284 population level of 2,000 and recorded the relative abundance of each of the two 

285 competing populations every 1,000 updates over the course of the 10,000 update trial. 
286 

287 In the second set of Phase 3 trials (Figure 1, Phase 3b), each of the starting ("pre-Phase 

288 2") and final ("post-Phase 2") Phase 2 prey populations were combined with a copy of 

289 the new Phase 3 predator population and reintroduced into a fresh environment. We then 

290 recorded the number of moves and turns executed and the usage of visual sensors ("look" 

291 instructions), as proxies for anti-predator responses. Each prey population was then 

292 introduced a second time in a separate evaluation in a predators absent PT environment. 

293 For each of these trials, we kept the number of prey at a constant 1,000 individuals via 

294 enforcement of population caps (i.e. a random prey was killed if a new birth would bring 

295 the prey population above 1,000) and minimum thresholds (below which predator attacks 

296 would fail). 
297 
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298 Statistical analysis 

299 To analyze the effects of SGV, EH, and PT outcomes of ecological competition (Phase 3; 

300 Figure 2), we used a linear mixed effects model (Appendix A) with SGV, EH, and PT as 

301 fixed effects, and a random intercept and slope across levels of all factors as random 

302 effects (Table S3), to account for non-independence of factor levels. For the response 

303 variable, a single observation was the mean relative abundance of a single treatment 

304 group in a given competition scenario (i.e. one of the two dots in each panel of Figure 2). 

305 All first-order interaction terms were included in the model to highlight specific trends. 

306 AIC indicated that a second-order interaction did not improve the model (AIC full model: 

307 -103.73; AIC reduced model: -1 12.87; Figure S2). 
308 

309 To generate response variables for models describing change in traits due to Phase 2 

310 evolution, we divided instruction counts by total instructions performed (e.g., number of 

311 moves / total lifetime instructions) to derive the proportion of instructions represented by 

312 each instruction type. Trait assays were conducted pre- and post-Phase 2. Thus, to 

313 determine the effect of evolution specifically during Phase 2, we subtracted pre-Phase 2 

314 instruction proportions from post-Phase 2 proportions, yielding the raw change in 

315 instruction proportions. All statistical analyses were performed on data scaled by the 

316 inverse of the variance to account for heteroscedasticity among treatment groups, though 

317 figures in the main text show unsealed data. 
318 

319 We used linear mixed effects models (Equation 1; Appendix B) fit by restricted 

320 maximum likelihood to determine the effects of SGV, EH, and Phase 3 predator 
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321 treatment (PT) on change in total instructions and in instruction proportions due to Phase 

322 2 evolution. SGV, EH, and PT were included as fixed effects. Since both predator 

323 treatment (PT) levels {predators present, predators absent) were applied to each 

324 population, we fitted replicate as a random intercept, and included a random slope across 

325 PT levels to account for non-independence of PT levels. Equation 1 describes the change 

326 in total prey instructions due to Phase 2 evolution. Models describing change in prey 

327 moves, turns, and looks varied in complexity, but were of the same general form 

328 (Appendix B). Models describing change in total instructions and changes in proportion 

329 moves, proportion turns, and proportion looks were of the same form, but varied in 

330 complexity. For each trait, we used a parametrically bootstrapped likelihood ratio test 

331 with for model selection, starting with the full model (Equation 1) and dropping 

332 interaction terms that did not improve the model fit (Appendix B, Table S4). 
333 

334 EQUATION 1 

335 y~fiov+ PiyPT +faPT +p 2 EH +frSGV +p 4 (PT *EH) +fr(PT xSGV) +fi 6 (EH x SGV) + 

336 p 7 (PT x EH x SGV) + e ijk , faj-Ntfoj, ofa, faij~N(faj, ojj, 
337 

338 where Poij are the random intercepts, and/?;/, are the randomly varying slopes across PT 

339 levels. We used a method recently developed by [29] to calculate marginal R 2 (proportion 

340 of variation accounted for by fixed effects alone) and conditional R 2 (proportion of 

341 variation accounted for by fixed and random effects together). R„ and for each model 

342 are given in Appendix B, Table S5. 
343 
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344 To test the effects of SGV and EH on predator attack rates during trait assays, we used a 

345 general linear model with SGV and EH included as fixed effects, and the raw difference 

346 in attack rates between pre- and post-Phase 2 evolution as the response variable. A single 

347 replicate population was excluded from analyses of change in total instructions, as this 

348 replicate had an unusually high pre -Phase 2 instruction count, was highly 

349 unrepresentative of prey populations in general, and greatly increased the variance of the 

350 clone SGV + no predator EH treatment group (Figure S7). 
351 

352 Software & Hardware 

353 We used Avida version 2.14 and the Heads-EX hardware [30] for all experiments. Avida 

354 did not originally support immortal, sterile predators as used in Phase 2 and Phase 3, so 

355 we added an additional option to Avida that, when set, causes predators to 'reset' (akin to 

356 being reborn) instead of reproducing or dying of old age. These modifications are 

357 available as a patch file against Avida 2.14, along with all analysis scripts, in the GitHub 

358 repository ( https :// github .com/ fishj ord/ avida_predation_scripts ) . 
359 

360 We conducted all experiments using the Michigan State University's High Performance 

361 Computing Cluster (HPCC). The MSU HPCC contains three general purpose compute 

362 clusters containing 2944 cores. Individual Avida runs were submitted as compute jobs to 

363 the HPCC general processing queue in parallel where possible. The Avida jobs took 

364 between 30 minutes and seven days to run, depending on the number of updates (2,000 to 

365 2,000,000 updates). 
366 
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367 We performed all statistical analyses and constructed all figures using the R statistical 

368 programming language version 3.0.2. Linear mixed effects models were constructed 

369 using the "lmer" function from the "lme4" package (last updated by [31]), and we used 

370 the "allEffects" function from the "effects" package (last updated by [32]) to extract 

371 marginalized fixed effects. To calculate marginal and conditional R 2 , we used the 

372 "r.squaredGLMM" function from the "MuMIn" package (last updated by [33]). We 

373 performed bootstrapped likelihood ratio tests using the "PBmodcomp" function from the 

374 "pbkrtest" package (last updated by [34]). 
375 

376 Results 

377 Ecological competitions 

378 Evolutionary history, i.e. historical exposure to predation in Phase 1, was the most 

379 important selective agent shaping competitive abilities of prey both in the presence and in 

380 the absence of a novel predator. Prey that evolved with predator EH were, in general, 

381 stronger competitors than those that evolved with no predator EH. In competitions 

382 between predator EH and no predator EH treatments, predator EH treatments always 

383 had the higher mean final relative abundance (Figure 2; Figure S2 D). Competitive 

384 exclusion of a predator EH lineage by a no predator EH lineage after 10,000 updates in 

385 competition was extremely rare (but see e.g. Figure 2, panel D3, in which 1 out of 30 

386 predator lineages was excluded). Competitive outcomes were not affected by the 

387 presence of a novel predator, regardless of SGV or EH (Figure 2; Figure S2 B, C). While 

388 there was no overall effect of SGV on mean final relative abundances, higher SGV seems 

389 to have conferred some benefit to no predator EH populations (Figure 2; Figure S2 B). 



17 



Downloaded from http://biorxiv.org/on September 18, 2014 

390 Similarly, there appears to be a subtle interaction between PT and EH, wherein, in the 

391 presence of a predator, the difference in relative abundances between predator EH and no 

392 predator EH populations during competition assays increases in favor of predator EH 

393 populations. 
394 

395 Evolution of prey traits during Phase 2 

396 Lifetime instruction counts for prey do not reflect higher rates of instruction execution, 

397 but rather longer lifespans of prey organisms. Changes in instruction counts varied 

398 greatly among treatments (PT x EH x SGV interaction; Figure S3). During Phase 2 

399 evolution, regardless of presence or absence of a novel predator during trait assays, 

400 predator EH populations generally evolved to increase total instructions executed per 

401 lifetime, except in the case of clone populations in the absence of predators (Figure 3). 

402 Conversely, no predator EH populations universally decreased in lifetime instruction 

403 counts. In predator EH populations, SGV had a markedly positive effect on the 

404 magnitude of the decrease in instruction counts, while SGV had little or no effect on 

405 instruction counts in predator EH populations. Prey generally changed more in total 

406 instructions executed with predators present than with predators absent (Figure 3). 
407 

408 Prey moves as proportions of total prey instructions universally decreased in predator EH 

409 populations as a result of Phase 2 evolution, and largely increased in no predator EH 

410 populations (EH effect, Figure 4; Figure S4). Effects of SGV on change in moves were 

411 subtle, though slightly less so with predators absent. There was no main effect of PT on 

412 change in moves (Figure 4; Figure S4). 
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413 

414 Prey turns as proportions of total prey instructions increased dramatically in predator EH 

415 populations, but changed very little in no predator EH populations (Figure 4; Figure S5). 

416 Effects of SGV were subtle, but were most pronounced in no predator EH populations: 

417 clone and intermediate SGV populations decreased slightly on average in proportion of 

418 moves, while high SGV populations increased (EH x SGV interaction, Figure S5). PT 

419 did not affect change in turns (no PT main effect, Figure 4; Figure S5), nor was there a 

420 main effect of SGV (Figure S5). 
421 

422 Proportions of instructions that were looks generally decreased in predator EH 

423 populations, but changed little in no predator EH populations; EH effects varied among 

424 levels of SGV, with the smallest EH effects in clone SGV populations (EH x SGV 

425 interaction, Figure 4; Figure S6). Main effects of SGV and PT were not statistically 

426 significant (Figure S6). 
427 

428 Effects of Phase 2 prey evolution on predator attack rates 

429 Phase 2 evolution led to a reduction in predator attacks in nearly all cases, and 

430 evolutionary history and standing genetic variation both affected predator attacks on prey 

431 (Figure 5; Table S2). High and intermediate SGV populations experienced qualitatively 

432 similar reductions in predator attacks as a result of Phase 2 evolution, and prey evolved 

433 with predator EH experienced similar reductions in predator attacks across SGV 

434 treatments (Figure 5). Only clone SGV + no predator EH populations were attacked more 

435 after Phase 2 evolution. 
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436 

437 Discussion 

438 Though the effects of evolutionary history (EH) and standing genetic variation (SGV) on 

439 the evolutionary potential of populations have been the subject of considerable conjecture 

440 [5,9], they have rarely been tested. Here, we were able to separately control for and vary 

441 the background levels of SGV and EH, and subsequently test the roles of each in shaping 

442 the evolvability of populations in a predator-prey context. 
443 

444 Role of evolutionary history 

445 Evolutionary history (the combined selective pressures under which populations evolve) 

446 likely determines the nature of mutations that fix in populations [9]. We found that 

447 having an evolutionary history with a predator improves a population's ability to evolve 

448 defenses in response to novel pressures from novel invading predator populations: in 

449 general, predator EH populations experienced greater changes in trait expression 

450 compared to no predator EH populations (Figure 4) when evolving defenses against 

451 invading (Phase 2) predators. This occurs regardless of the level of SGV, suggesting that 

452 the importance of EH eclipses that of evolvability originating from SGV. Furthermore, in 

453 pairwise competitions (Figure 2), populations with predator EH demonstrated 

454 competitive superiority over no predator EH populations, similar to findings of [18]. If a 

455 novel predator was present during the competition, the advantage was even more 

456 pronounced (also see Figure S2 C), with the vast majority of predator EH populations 

457 competitively excluding their no predator EH competitors (Figure 2). Additionally, when 

458 introduced into separate ecological trials with the novel predator, populations that 



20 



Downloaded from http://biorxiv.org/on September 18, 2014 

459 evolved with predator EH suffered fewer attacks than populations evolved with no 

460 predator EH (Figure 5). 
461 

462 Because both types of EH populations were further evolved in the presence of predators 

463 (in Phase 2) before they were competed against each other, the ecological advantage 

464 enjoyed by predator EH populations in competition with no predator EH populations 

465 was due, in part, to their historic evolution of anti-predator traits (i.e. in Phase 1), even 

466 though that initial evolution was in response to an entirely different predator population. 
467 

468 Both the mutation rates used for the evolution of phase 1 predator EH populations and the 

469 duration of the evolutionary trials will have influenced the extent of evolutionary history 

470 in each trial (i.e. the degree to which prey populations were adapted to predators). 

471 Because we did not vary Phase 1 mutation rates or duration, we cannot comment on how 

472 experimentally altering the strength of evolutionary history could impact the observed 

473 trends. 
474 

475 Role of standing genetic variation 

476 Many studies have pointed to SGV as a key component of a population's evolutionary 

477 potential [6]. When there is greater genetic variation in a population, there is more raw 

478 material upon which natural selection can act. Thus, we expected SGV alone to play a 

479 significant role in the rapid evolution of anti-predator traits. However, we found only 

480 small effects of historic SGV in determining the final outcomes of ecological 

481 competitions (Figure 2). Furthermore, the effects of SGV on the evolution of specific 
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482 anti-predator traits (moves, turns, looks) were detectable only within EH treatments 

483 (Figure 4). 
484 

485 Standing genetic variation (SGV) vs Evolutionary history (EH) 

486 Across the predator EH populations, we found no detectable effects of SGV on anti- 

487 predator trait values (Figure 4) or competitive outcomes (Figure 2). We suggest that a 

488 lack of an SGV effect on predator EH populations may be due, in part, to the historical 

489 predator weeding out unfit prey genotypes. E.g., predators reduced prey diversity (Figure 

490 SI; Table SI) because they operate as a very strong agent of selection against less fit 

491 phenotypes. However, we would expect SGV to have an effect within no predator EH 

492 populations (as seen in Figure 4). Given the lack of an evolutionary history with 

493 predators, in no predator EH populations, the odds of one sampled genotype {clone SGV) 

494 being able to quickly acquire a beneficial mutation are low. On the other hand, if the full 

495 suite of all discovered genotypes {high SGV) remains in the population, there is a greater 

496 potential for providing precursory material for rapidly generating anti-predator traits. 

497 However, if EH leads to a large number of genotypes having the prerequisites for rapid 

498 adaptive response, the odds of sampling a beneficial genotype are high, and differences in 

499 SGV will not significantly affect evolvability. 
500 

501 Indeed, as we have shown, in order to adaptively address novel predation pressures, it 

502 appears to be easier to modify historically realized and relevant adaptations than to 

503 repurpose unrelated genes, or evolve effective traits de novo. Thus, the greater 

504 evolvability of the predator EH populations would have arisen out of their ability to 
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505 adapt anti-predator strategies across predator contexts. More broadly, evolvability may 

506 commonly be determined by a population's ability to use strategies across environmental 

507 contexts, i.e. their adaptive complexity [35]. 
508 
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596 Figure legends 

597 Figure 1 

598 Design of the experiment to test the effects of SGV and EH on evolution of anti-predator 

599 traits. A) 9 identical prey are injected to initialize two sets of 30 base populations. B) 

600 Over the course of 2 million updates of evolution, populations diversified, including the 

601 de novo evolution of predators in half of these EH populations (top). C) One fully 

602 evolved prey population from each EH treatment (with and without evolving predators) 

603 was used to create 30 identical sets of high, intermediate, and no diversity (clone) SGV 

604 populations (totaling 2 X 30 X 3 = 180 populations), with clones of the best evolved 

605 predator from Phase 1 added to each. D) After a further 2 million updates of evolution at 

606 a low mutation rate, replicate populations converged on an intermediate level of diversity, 
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607 but SGV and EH differences result in trait variation among populations. E) To evaluate 

608 the effectiveness of evolved anti-predator traits, each of the 180 fully evolved populations 

609 from D were introduced, in turn, into separate competition assays with one randomly 

610 selected replicate from each of the other EH X SGV combination populations, both in the 

611 presence and in the absence of a new novel predator. F) Separately, prey rates of 

612 executing "anti-predator" actions (moves, turns, looks) were recorded in the presence and 

613 absence of predators, to quantify expressed behavioral responses to predators for all fully 

614 evolved populations . 
615 

616 Figure 2 

617 Randomly paired populations in competitions for all SGV and EH treatment 

618 combinations in the absence (gray boxes, above diagonal) and presence (black boxes, 

619 below diagonal) of a novel predator (Phase 3). EH had a far greater effect on competitive 

620 outcomes than SGV in both the presence and absence of predators, with predator EH 

621 populations competitively dominant in most cases. Y-axis: relative abundance of each 

622 lineage in competition. Colors for each lineage indicate SGV and EH treatment history. 

623 Black: clone, predator, Gray: clone, no predator, Blue: intermediate, predator, Red: 

624 intermediate, no predator, Gold: high, predator, Pink: high, no predator. Dots indicate 

625 mean relative abundance for lineages of each treatment after 10,000 updates of 

626 competition. 
627 

628 Figure 3 
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629 Change in total prey instructions executed (start of Phase 2 minus end of Phase 2), 

630 reflecting change in gestation time (or lifespan) resulting from Phase 2 evolution. 

631 "Predator present/ absent" refers to the presence or absence of a novel predator during 

632 trait assays before and after Phase 2 evolution. Total instructions either increased or 

633 decreased very little in predator EH populations, and universally decreased in no 

634 predator EH populations. Bars are ± 95% CI. 
635 

636 Figure 4 

637 Change in (top) prey moves, (middle) prey turns, and (bottom) prey looks as proportions 

638 of total prey instructions, resulting from Phase 2 evolution with a novel predator. 

639 "Predator present/absent" refers to the presence/absence of the Phase 3 novel predator 

640 during the trait assays. Traits in predator EH populations generally evolved to a greater 

641 extent than did traits in no predator EH populations, and traits often evolved in different 

642 directions between the two EH treatments. Bars are ± 95% CI. 
643 

644 Figure 5 

645 Change in attacks in behavioral assays pairing the Phase 3 predator with pre- and post- 
646 Phase 2 prey populations. Predator attack frequencies on predator EH populations 

647 universally decreased, and generally changed more than did attack frequencies on no 

648 predator EH populations. Bars are ± 95% CI. 
649 

650 Supporting figure legends 

651 

652 Figure SI 
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653 Shannon's diversity across 200,000 updates of Phase 2 evolution in (A) high, (B) 

654 intermediate, and (C) clone SGV treatments. Filled circles represent no predator EH 

655 treatments, open circles represent predator EH treatments. Bars are ± 95% CI. 
656 

657 Figure S2 

658 Marginalized fixed effects of a generalized linear mixed effects model describing effects 

659 of SGV, EH, and PT on mean relative abundance of a single treatment group from each 

660 competition scenario depicted in Figure 2. A. SGV x PT interaction plot; B. SGV x EH 

661 interaction plot; C. EH x PT interaction plot. Bars are ± 95% CI. EH was by far the most 

662 important factor explaining competitive outcomes; predator EH groups were universally 

663 superior to no predator EH groups in ecological competition. Formal model, random 

664 effects output, and marginal and conditional R 2 can be found in Appendix A. 
665 

666 Figure S3 

667 Marginalized fixed effects of a generalized linear mixed effects model of scaled change 

668 in total prey instructions across Phase 2 evolution as a function of SGV (A. clone; B. 

669 intermediate; C. high), EH, and PT. Bars are ± 95% CI. EH had the greatest effect on 

670 change in prey instructions, though PT also affected change in instructions, especially in 

671 clone SGV populations. Formal model, random effects output, and marginal and 

672 conditional R 2 can be found in Appendix B. 
673 

674 Figure S4 
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675 Marginalized fixed effects of a generalized linear mixed effects model of scaled change 

676 across Phase 2 evolution in the proportion of prey instructions that were moves as a 

677 function of A. SGV and PT and B. EH. Bars are ± 95% CI. Both EH and SGV affected 

678 change in proportion moves, with SGV increasing more dramatically across SGV 

679 treatments with predator present. Formal model, random effects output, and marginal and 

680 conditional R 2 can be found in Appendix B. 
681 

682 Figure S5 

683 Marginalized fixed effects of a generalized linear mixed effects model of scaled change 

684 across Phase 2 evolution in the proportion of prey instructions that were turns as a 

685 function of A. PT and B. SGV x EH. Bars are ± 95% CI. Predator EH populations 

686 increased turns as a results of Phase 2 evolution; change in turns was marginal for no 

687 predator EH populations, and SGV affected predator and no predator populations 

688 differently. PT had no detectable effect. Formal model, random effects output, and 

689 marginal and conditional R 2 can be found in Appendix B. 
690 

691 Figure S6 

692 Marginalized fixed effects of a generalized linear mixed effects model of scaled change 

693 across Phase 2 evolution in the proportion of prey instructions that were looks as a 

694 function of A. PT and B. SGV x EH. Bars are ± 95% CI. Clone and high SGV 

695 populations with predator EH increased looks as a result of Phase 2 evolution and change 

696 in turns was marginal to negative for no predator EH populations. PT had no detectable 
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697 effect. Formal model, random effects output, and marginal and conditional R" can be 

698 found in Appendix B. 
699 

700 Figure S7 

701 Total instructions executed by each replicate prey population before and after Phase 2 

702 evolution. The slope of each line represents the evolutionary change in instructions 

703 executed per lifetime (i.e. gestation time). The same replicate populations were assayed 

704 in the absence (top) and in the presence (bottom) of the Phase 3 predator. Arrows indicate 

705 the replicate population that was excluded from all statistical analyses, as it is not 

706 representative of its treatment group {clone SGV, no predator EH). 
707 

708 
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1 Supporting Information (SI) 

2 

3 Genetic diversification during Phase 2 evolution 

4 Over the 200,000 updates of Phase 2 evolution, clone SGV populations increased in 

5 diversity (Shannon's Diversity Index), and intermediate and high SGV populations 

6 decreased in diversity, with all trends stabilizing near the 40,000 th update (Figure SI). 

7 Both standing genetic variation (SGV) and evolutionary history (EH) affected diversity 

8 values at the end of Phase 2 evolution. Diversity was similar across SGV treatments for 

9 predator EH populations, but was higher in intermediate and high than in low SGV 

10 populations for no predator EH populations. Shannon's Diversity was higher over all for 

11 no predator than for predator EH populations (Figure SI; Table SI). 
12 

13 Supplemental tables (S2) 
14 

15 Table SI Analysis of Variance (calculated using the anova() function in the R 3.0.2 base package) on a 

16 linear model*' testing the effects of SGV and EH and their interaction on Shannon's diversity among 

1 7 populations at the 200,000 th update of Phase 2 evolution. 



Effect 


DF 


Mean squares 


F-value 


P-value 


SGV 


2 


1.45 


2.45 


0.09 


EH 


1 


34.74 


59.72 


< 0.0001 


SGV x EH 


2 


2.65 


2.28 


0.11 


Error 


174 


0.582 







18 *Multiple R 2 = 0.45 

1 9 fModel: y ~ p 0 + faSG V + p 2 HC + p 3 (SG V x HQ + e 
20 

21 Table S2 Analysis of Variance (calculated using the anova() function in the R 3.0.2 base package) on a 

22 linear model* ' testing the effects of SGV and EH and their interaction on change in predator attacks as a 

23 result of Phase 2 evolution. 
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Fixed effect 


DF 


Mean squares 


F-value 


P-value 


SGV 


2 


1.60 x 10 15 


19.89 


< 0.0001 


EH 


1 


4.30 x 10 15 


53.51 


< 0.0001 


SGV x EH 


2 


8.41 x 10 13 


1.05 


0.35 


Error 


354 


8.04 x 10 13 







24 *MultipleR 2 =0.20 

2 5 fModel: y ~ p„ + fi,SG V + fi 2 EH + p 3 (SG V * EH) + e 
26 

27 Appendix A (S3) 

28 

29 Linear mixed effects model describing effects of SGV, EH, and PT on mean final relative 

30 abundances of treatment groups in pairwise ecological competitions (Figure 2; Figure 

31 S2). 
32 

33 Model A. 1: 

34 y ~ fioff + PnjPT + p 2ij EH + jS 3ij SGV +JS 4 PT + fi 5 EH + fcSGV + JS 7 (SGV x EH) + 

35 fa(SGV x PT) + fi g (EH x PT) + s ijh 

36 Poij-Ntfoj, o} Q ), Pnj-NfPu, a} J, Pnj-N&j, oj 2 ), fi 3v ~ N0 3J , tf 3 ) 
37 

38 Poij are the randomly varying intercepts, /3nj are the randomly varying slopes across PT 

39 levels, /?2iy are the randomly varying slopes across EH levels, and are the randomly 

40 varying slopes across SGV levels. 
41 

42 Table S3 Random intercepts and slopes (linear contrasts) fit in Model A. 1 with parametric bootstrap 95% 

43 confidence intervals.* Different replicate ID's were assigned by treatment combination for each factor. 

44 =0.51; R 2 C = 1.0. 



Factor 


Coefficient 


Std. dev. 


Lower 95% 


Upper 95% 


SGV 




slope (intermediate / clone) 


0.030 


0.030 


0.030 



2 
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slope (high / clone) 


0.030 


0.030 


0.030 




intercept 


0.18 


0.17 


0.18 


EH 




slope 


0.038 


0.038 


0.039 




intercept 


0.031 


0.030 


0.031 


PT 




slope 


0.015 


0.015 


0.015 




intercept 


0.0088 


0.0087 


0.0089 


Residuals 




0.0018 


0.0017 


0.0018 



45 *Some confidence interval values are the same as standard deviations, due to rounding. 
46 

47 Appendix B (S4) 

48 

49 Linear mixed effects models describing change in prey instructions across Phase 2, 

50 compared using LRT with a parametric bootstrap. Models are nested, with Model 1 being 

5 1 the full model (all first and second-order interactions included), and Model 6 being the 

52 simplest (no interactions included). Models 2-5 include some (Models 3-5) or all (Model 

53 2) first-order interactions. Interaction terms included were those determined to be 

54 statistically significant (parametric bootstrap confidence intervals around coefficients did 

55 not overlap zero). LRT outputs, AIC for each model, and model selection processes are 

56 given in Table S3. 
57 

58 Model B.l (full model): 

59 y ~ fa + fiiijPT +fiiPT +p 2 EH +0 3 SGV +p 4 (PT xEH) +p 5 (PT xSGV) + 

60 p 6 (EH x SGV) + p 7 (PT xEHx SGV) + e iJh 

61 Poij~N(p 0] , tfa), Puj-NfPu, aj& 
62 

63 Model B.2: 

64 y ~ Poij + PnjPT +PiPT +p 2 EH +p 3 SGV +p 4 (PT *EH) +p 5 (PT *SGV) + 
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65 p 6 (EH x SGV) + s ijh 

66 Poij ~ N(p 0j , tfo), Phj ~ 
67 

68 Model B.3: 

69 y~Poij + PiijPT + PiPT + p 2 EH + p 3 SG V + p 4 (EH x SGV) + s ijh 

70 Poij ~ iV(% o^, Phj ~ iV(% o-|v> 
71 

72 Model B.4: 

73 j ~ Poij + PnjPT +piPT + p 2 EH + p 3 SGV + p 4 (PT x SGF) + 

74 ^ ~ N(P 0] , tfo), Phj ~ N(Pu, ofa 
75 

76 Model B.5 (fully reduced model): 

77 y ~ Poij + PnjPT +P,PT + p 2 EH + p 3 SGV + e ijh 

78 ^ ~ A^, o-|o>, ^ ~ 
79 

80 Pay are the randomly varying intercepts, and P nj are the randomly varying slopes across 

81 PT levels. 
82 

83 Table S4 AIC, LRT output, and model selection processes for linear mixed-effects models describing 

84 change in A. total prey instructions; B. proportion moves; C. proportion turns and D. proportion looks. 



Compared 


AIC 


X 2 


df 


LRT P-val 


PBtest P-val 


Mod. preferred 


A. Total instructions 


B.l &B.2 


678.09 


692.23 


18.14 


2 


0.0001 


0.001 


1, 1 chosen 


B. Proportion moves 


B.l &B.2 


370.83 


370.82 


3.98 


2 


0.14 


0.15 


none, 2 chosen 


B.2 & B.4 


370.82 


367.76 


1.99 


3 


0.57 


0.55 


none, 4 chosen* 


C. Proportion turns 
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85 
86 
87 
88 
89 
90 



B.l & B.2 


587.37 


586.37 


3.98 


2 


0.14 


0.16 


none, 2 chosen 


B.2 & B.4 


586.37 


589.67 


9.30 


3 


0.03 


0.03 


2, 2 chosen 


D. Proportion looks 


B.l & B.2 


598.84 


595.98 


1.14 


2 


0.57 


0.59 


none, 2 chosen 


B.2 & B.3 


595.98 


592.18 


2.20 


3 


0.53 


0.53 


none, 3 chosen 


B.3 & B.5 


592.18 


594.18 


6.00 


2 


0.05 


0.05 


3, 3 chosen 



*A11 remaining first-order interactions significant, so no further model comparison. 

Table S5 Random intercepts and slopes (linear contrasts) between levels of PT with parametric bootstrap 
95% confidence intervals for models describing change in A. total prey instructions; B. proportion moves; 
C. proportion turns; and D. proportion looks. Also included are marginal and conditional R 2 and residual 
standard deviations with parametric bootstrap 95% CI. 



Trait 


Coefficient 


Std. dev. 


Lower 95% 


Upper 95% 


A. Total instructions: RL = 0.57; Ri = 0.82 




slope 


0.46 


0.38 


0.52 




intercept 


0.29 


0.21 


0.35 




residuals 


0.39 


0.34 


0.41 


B. Proportion moves: Ri, = 0.56; Ri = 0.97 




slope 


0.61 


0.54 


0.67 




intercept 


0.02 


0.004 


0.06 




residuals 


0.17 


0.15 


0.19 


C. Proportion turns: R^ = 0.23; R^= 1.0 




slope 


0.81 


0.72 


0.88 




intercept 


0.37 


0.32 


0.40 




residuals 


0.00044 


0.00040 


0.00048 


D. Proportion looks: Ri, = 0.10; R 2 C = 1.0 




slope 


0.96 


0.83 


1.03 




intercept 


0.27 


0.23 


0.31 




residuals 


0.13 


0.12 


0.14 
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